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Abstract. We show that in the system of two fermions interacting by scalar exchange, the solutions 
for J'^=0+ bound states are stable without any cutoff regularization, for values of the coupling 
constant a below a critical value a^ . This latter is calculated from an eigenvalue equation. 



^ : INTRODUCTION 

m 

CN . The Yukawa model, describing a system of two fermions interacting by scalar exchange 

^j^int _ gi|rv|/(|) ('')), is instructive for studying the relativistic bound states and for de- 
veloping the renormalization methods. It is also a main ingredient in building the A^A^ 
interaction, which contains an important contribution of scalar meson exchanges. The 
bound state problem and the renormalization in the Yukawa model were studied [1] in 
the framework of standard light-front dynamics [2]. The relativistic two-nucleon wave 
Q^' functions have been also calculated perturbatively [3] in the explicitly covariant ver- 

^ . sion of light-front dynamics [4], where the state vector is defined on the invariant plane 

Qdx = with (0^ = (see for review [5]). In this work, the Bonn A^A^ model was used 
with the corresponding form factors [6] and the problem of cut-off dependence was not 
analyzed. 

In reference [7], we investigated the stability of the bound states relative to the high 
momentum contributions of the kernel, when the cutoff tends to infinity. Below we 
present the results of our study and compare them with those obtained in [1]. 



THE CUTOFF DEPENDENCE OF THE BINDING ENERGY 

We consider the two fermion wave function with total angular momentum 7 = 0+. In 
the fermion spin indices, it is a 2 x 2 matrix determined, due to the parity conservation, 
by two independent elements - the spin components fi=\,2- Since the wave function 
is defined on the light- front plane CO-x = 0, components fi{k,Q) depend not only on 
the relative momentum k, but also on the angle between k and h = (b/\(b\. The 
system of equations for /, contains a 2 x 2 matrix kernel Kij, calculated by using 
the explicitly covariant light-front graph techniques [4, 5]. These equations and the 
analytical expressions of kernels are given in [7]. Though their form differs from the 
ones used in [1], we have shown that they are strictly equivalent. 



Let us consider the equations on the finite interval <k < kmax- The dependence of 
their solutions on the cutoff kmax in the limit kmax ^ 0° is determined by the kernels 
asymptotics. The kernel K22 is repulsive and cannot generate a collapse, whereas ^n 
is attractive. Therefore, to investigate the stability, we can consider the one channel 
problem for the component /i, which satisfies the equation: 

[M^-4ik^ + m^)]Mk,Q) = ^jKnik,Q;k',Q')Mk'M)^, (D 

where Ej^i = y/nf- + k'^. 

Our analysis uses the fact that atk^ 0° the integral in the r.h.s. of (1) is dominated by 
the region k' oc k, i.e. k' ^ °° with a fixed ratio k' /k = Y [8]. One can therefore replace in 
(1) both wave function and kernel by their asymptotics, which have the form [7] : 

/■rt.^-M£l ir __^f Vy^n{z,z',y), ifY<l .^. 

^H'^'^^-^2+p' ^11- m2\ Aii(z,z',l/Y)/V?, ifY>l ^^ 

with < (3 < 1 and 



A ( ' )= f^'' "^^ 2Y(l-zzO-(l+f)v^T^?v^T^Fcos(^ 

where a = g^/4-n and z = cosG. By setting Aii(z,z',y) = 1 and a = 2nma' in (2), the 
asymptotics of ^n becomes identical to asymptotics of the momentum space kernel 
corresponding to the non-relativistic potential —a'/r^. As it is well known [9], there 
exists for this potential a critical coupling constant a'. = l/(4m), that corresponds to 
Uc = Tl/l. The inspection of (3) shows that in the Yukawa model the function An is 
smaller than for the —<x'/r^ potential: < Aii(z,z',y) < 1. Therefore in this model, 
one can expect a larger critical coupling constant i.e. QLc > ll/2, what is confirmed by 
numerical calculations. 

Substituting (2) into equation (1), we obtain for ^1 (z) [8]: 

j_ dz'H^{z,z)hi{z')=Xhi{z) (4) 

with ?i = 1/a and 

//p(z,z') = /^'^^iifc2',Y) cosh(|3logY). (5) 

Equation (4) is an eigenvalue equation for X, parametrized by (3. It provides the relation 
between the coupling constant a and (3, determining the power law (2) of the wave 
function asymptotics. The r.h.s. of equation (1) becomes divergent for |3 < 0. Hence, the 
equation |3(at) = determines the critical coupling constant ttc. 



RESULTS 



In the numerical calculations, the constituent masses were taken equal to m=l and 
the mass of the exchanged scalar is /j=0.25. The numerical solution cx((3), found from 
equation (4) with the function Aii(y,z,z') given by (3), is plotted in figure 1. The critical 
coupling constant is obtained for (3 = for which the eigenvalue is Xc = 0.269. It 
corresponds to a^ = Ifkc = 3.72, as shown in figure 1 at (3 = 0. 
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FIGURE 1. Function a(p) determined by eq. (4). The critical coupling constant is: a^ = a(P = 0) = 
3.72. The values discussed in the text: a = 1.096, p == 0.819 and a = 2.480, p = 0.548 are on the curve. 

In figure 2, we have plotted the mass square M^ of the two fermion system, found 
from (1), as a function of the cutoff kmax for two fixed values of the coupling constant 
below and above the critical value a^ = 3.72. One can see two dramatically different 
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FIGURE 2. Cutoff dependence of the binding energy in the J — 0^ state, in the one-channel problem 
(/i), for two fixed values of the coupling constant below and above the critical value. 

behaviors depending on the value of the coupling constant a. For a = 3, i.e. a < ttc, the 
result is convergent. On the contrary, for a = 4, i.e. a > ttc, the result is clearly divergent: 
M^ decreases logarithmically as a function of kmax and becomes even negative. Though 



the negative values of M^ are physically meaningless, they are formally allowed by 
equations. The first degree of M does not enter neither in the equation nor in the kernel, 
and M^ crosses zero without any singularity. 

We have examined the asymptotical behavior of the wave function /i {k, z) and found 
that it very accurately follows the power law (2) with the power |3(a) given in figure 1. 
For instance for a binding energy B = 2m — M = 0.05 (a = 1 .096) a direct measurement 
in the numerical solution plotted in figure 3 gives |3 = 0.820 ±0.002 whereas the solution 
of equation (4) for the corresponding a gives |3 = 0.819. The same kind of agreement 
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FIGURE 3. Asymptotic behavior of the 7 = 0+ wave function components /i,2 for B=0.05, a=1.096, 
/j=0.25. The slope coefficients are pi — 0.820 and P2 ~ 0. 

was found for B = 0.5 (a = 2.480, (3 = 0.548 ± 0.002). 

We conclude that the 7 = - or (1+, 2—) state in the classification [1] - is stable (i.e. 
convergent relative to the cutoff kmax -^ °°) for coupling constant a below the critical 
value at = 3.72. In this point, our conclusion differs from the one settled in [1], where 
it was stated that the integrals in the r.h.s. of the equations diverge logarithmically with 
cutoff. Above the critical value the integrals indeed diverge and the system collapses. 

In the 7=1 state the system is found to be always unstable, as pointed out in [1]. 

Thus, by an analytical method, confirmed by numerical calculations, we have shown 
that the Yukawa model is not cutoff dependent for coupling constant below a critical 
value. The results obtained should be taken into account for instance in the renormaliza- 
tion procedures. 
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